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We create and study persistent currents in a toroidal two-component Bose gas, consisting of 87 Rb atoms in 
two different spin states. For a large spin-population imbalance we observe supercurrents persisting for over two 
minutes. However we find that the supercurrent is unstable for spin polarisation below a well defined critical 
value. We also investigate the role of phase coherence between the two spin components and show that only the 
magnitude of the spin-polarisation vector, rather than its orientation in spin space, is relevant for supercurrent 
stability. 
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Persistent currents are a hallmark of superfluidity and su- 
perconductivity, and have been studied in liquid helium and 
solid state systems for decades. More recently, it became pos- 
sible to trap an atomic Bose-Einstein condensate (BEC) in a 
ring geometry |Q}[8) and induce rotational superflow in this 
system O O [HI - This offers new possibilities for fundamen- 
tal studies of superfluidity in a flexible experimental setting. 
Both long-lived superflow (HO and quantised phase slips cor- 
responding to singly-charged vortices crossing the superfluid 
annulus have been observed O [8) . 

So far experiments on persistent currents in atomic BECs 
were limited to spinless, single-component condensates. Ex- 
tending such studies to multi-component systems, in particu- 
lar those involving two or more spin states I9UTT1. is essential 
for understanding superfluids with a vectorial order parame- 
ter and for applications in atom interferometry [Q21 03) • Per- 
sistent flow in a two-component Bose gas has been studied 
theoretically Q3HI3 but many issues remain open. Even the 
central question of whether, and under what conditions, this 
system supports persistent currents has not been settled. 

In this Letter, we study the stability of supercurrents in a 
toroidal two-component gas consisting of 87 Rb atoms in two 
different spin states. For a large spin-population imbalance we 
observe superflow persisting for over two minutes and limited 
only by the atom-number decay. However at a small popula- 
tion imbalance the onset of supercurrent decay occurs within a 
few seconds. We demonstrate the existence of a well defined 
critical spin-polarisation separating the stable- and unstable- 
current regimes. We also study the connection between spin 
coherence and superflow stability, and show that in our sys- 
tem only the modulus of the spin-polarisation vector is rel- 
evant for the stability of the supercurrent. The existence of 
a critical population imbalance was anticipated in Refs. [TT3T— 
[TTl . but quantitative comparison with our measurements will 
require further theoretical work. 

Our setup is outlined in Fig. [TJa). We load a BEC of 
N w 10 5 atoms into an optical ring trap of radius 12 /mi, 
created by intersecting a 1070 nm "sheet" laser beam and an 
805 nm "tube" beam (7). The sheet beam confines the atoms 
to the horizontal plane with a trapping frequency of 350 Hz. 
In-plane, the tube beam confines the atoms to the ring with a 
trapping frequency of 50 Hz. The trap depth is about twice 



the BEC chemical potential, no/h « 0.6 kHz, and varies az- 
imuthally by < 10%. 

Our tube trapping beam is a Laguerre- Gauss LG 3 laser 
mode in which each photon carries orbital angular momentum 
3h. We use the same beam to induce a supercurrent via a two- 
photon Raman process GJ [71 H8) . We briefly (200 /is) pulse 
on an auxiliary TEM o Gaussian beam, copropagating with 
the LG beam, to transfer all atoms between two spin states, |t) 
and \l) [Fig.[ljb)]. Each atom absorbs angular momentum 3h 
from the LG beam and we thus create a (single-component) 
current corresponding to a vortex of charge q = 3 trapped at 
the ring centre. Such current can persist for over a minute and 




FIG. 1 : (color online) Preparation and detection of supercurrent in 
a two-component gas. (a) The ring trap is formed by a horizontal 
"sheet" beam and a vertical Laguerre-Gauss (LG) "tube" beam. B is 
the external magnetic field, (b) Supercurrent is induced by a Raman 
transfer of atoms between two spin states, |t) and ||), using the LG 
beam and an auxiliary Gaussian (G) beam. During the transfer each 
atom absorbs 3h of angular momentum from the LG beam. Two- 
component gas is created by coupling |f ) and |^) states with an RF 
field. The characteristic rotational energy is E r /h « 0.4 Hz. (c) 
Time-of-flight image of the atoms, with spin states separated using 
a Stern-Gerlach gradient. The rotational state q is deduced from the 
radius R characterising the central hole in the density distribution. 
The image shown was taken after t = 4 s of rotation; the longitudinal 
spin-polarisation is P z = 0.44 and q = 3 for both spin states, (d) 
Histogram of « 900 measurements of R at various P z and t. 
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FIG. 2: (color online) Single- vs. two-component supercurrent. (a) 
In a pure |f) state (P z = 1) supercurrent persists for over two min- 
utes, with no phase slips occurring for « 90 s. (b) At P z = the 
first phase slip occurs within 5 s and we observe no rotation beyond 
20 s. (c) Total atom number decay for P z = 1 (open symbols) and 
P z = (solid symbols). Dashed lines are double-exponential fits. 



decays in quantised q — » q — 1 steps, corresponding to 27r 
phase slips in the BEC wave function Q. 

The |t) and \\) states also define the spin space for our two- 
component experiments. To create a two-component current 
we prepare a pure \q = 3,|) state and then couple |t) and 
||) by a radio-frequency (RF) field, which carries no orbital 
angular momentum and does not affect the motional state of 
the atoms. The |t) and \\) are two F = 1 hyperfine ground 
states, mp = 1 and 0, respectively. The = -1 state 
is detuned from Raman and RF resonances by the quadratic 
Zeeman shift in an external magnetic field B of 10 gauss. 

After preparing a rotating (q = 3) cloud in a specific spin 
state, we let it evolve in the ring trap for a time t and then 
probe it by absorption imaging after 29 ms of time-of-flight 
expansion. We separate the two spin components with a Stern- 
Gerlach gradient and directly measure the longitudinal spin- 
polarisation P z (7V t - A^)/(A^ + N±), where 7V t (N±) 
is the number of atoms in the |t) (||)) state [Fig.[ljc)]. The 
rotational state, < q < 3, is seen in the size R of the central 
hole in the atomic distribution [7 ], arising due to a centrifugal 
barrier 0. As shown in Fig. [TJd), the R values are clearly 
quantised 00, allowing us to determine q with > 99% fi- 
delity E2. 

In Fig. [2] we illustrate the dramatic difference between su- 
perflow stability in a P z = 1 single-component gas and a 
P z = two-component system. The two different P z states 
are created, respectively, by a (140 /is) tt and a (70 /is) 7r/2 
RF pulse at t = 0. In the pure |t) state [Fig.[2ja)] the current 
persists for over two minutes, with the BEC always remaining 
in the q = 3 state for « 90 s. In contrast, at P z = [Fig.[2jb)] 
the first phase slip occurs within 5 s and the current completely 
decays within 20 s. During the decay we always observe the 
two spin components to be in the same q state. 

Supercurrent stability generally depends on the number of 
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FIG. 3: (color online) Supercurrent stability in a partially spin- 
polarised gas. The statistically averaged supercurrent state, (q), of 
the majority spin component is shown as a function of P z and the 
evolution time t. The contour plot is based on w 1600 measure- 
ments of q(P Z jt). The transition between stable- and unstable- 
current regimes occurs at 0.6 < P z < 0.7. In the stable regime 
the current eventually decays due to the atom-number decay. 

condensed atoms 13 and at P z =0 the atom number per 
spin state is halved. However, from the TV-decay curves in 
Fig. [2jc) we see that this alone cannot explain the difference 
in superflow stability. At P z = 1 rotation still persists for 
N « 10 4 while at P z = it stops already at N > 4 x 10 4 . 
Moreover, if we apply a tt/2 RF pulse at t = but then imme- 
diately remove all the |t) atoms from the trap with a resonant 
light pulse, the current again persists for over a minute. This 
unambiguously confirms that in Fig. |2jb) the superflow is in- 
hibited by the presence of both spin components. 

We now turn to a quantitative study of the supercurrent sta- 
bility as a function of the spin-population imbalance (Fig. [3]). 
We tune P z by varying the length At of the RF pulse applied 
at t = 0, and measure the q state of the majority (|t)) spin 
component as a function of t. Whenever the radius R is fit- 
table for the minority component we get the same q for both 
spin components in > 99 % of cases. However for N± < 10 4 
we cannot determine q for the minority component. 

Based on « 1600 measurements of q(P z ,t), in Fig. [3] we 
reconstruct the complete current stability diagram for < 
P z < 1 l20l . The contour plot of (q(P z ,i)) is obtained by 
spline interpolation through a 3D mesh of data points with in- 
teger q values. The blue-shaded region corresponds to rotation 
times for which no phase slips occur. We clearly distinguish 
two qualitatively different regimes. For large P z the superflow 
is fundamentally stable and limited only by the atom-number 
decay; for low P z the current starts to decay within a few sec- 
onds. A sharp transition between the two regimes occurs at 
0.6 <P Z < 0.7. 

To fully understand these observations, we need to distin- 
guish a coherent superposition of |t) and \\) states from an 
incoherent mixture. The RF pulse at t = corresponds to ro- 
tation around the y axis on a Bloch sphere and puts the BEC in 
a superposition state \0) = sin(#/2) |t) + cos(#/2) ||). Here 
= flRAt, where is the RF Rabi frequency. In this state 
P z = — cos(#) but the gas is still fully spin-polarised; the po- 



3 




t=60s 



FIG. 4: (color online) Role of spin-coherence in superflow stability, 
(a) An RF pulse at t = creates a spin- superposition state in which 
P z < 1 but P — 1. As the superposition decoheres, transverse 
spin-polarisation decays and P — >• P z . (b) Transverse-polarisation 
decay (see text). Double-exponential fit (solid red line) gives the 
decay function /(£). (c) Adiabatic dressing of the spin state. In 
presence of a resonant RF field 5 e ff oc y and the dressed P z = 
state \y) is stable against decoherence. (d) In the dressed \y) state 
the supercurrent is also stable. 

larisation vector is P = (sin(0), 0, — cos(0)) and P = |P| = 
1. Subsequently the spin superposition decoheres, due to both 
intrinsic spin diffusion [21 ] and small magnetic field inhomo- 
geneities [22 ]. P z is a constant of motion but the transverse 
polarisation decays and P — > P z [Fig.[4ja)]. 

We study the transverse-polarisation decay in a Ramsey- 
type experiment. Starting in the ||) state we apply two tt/2 
RF pulses separated by time t and then measure P z . The first 
pulse creates a purely-transverse P = (1,0,0) and the second 
one maps the decaying P into P z l23ll . As seen in Fig.|4|b), 
we observe a very long spin-coherence time, £ co h ~ 10 s. This 
means that in the unstable regime in Fig. [3] phase slips occur 
already at t < £ co h, when we cannot equate P and P 2 . 

We also perform a complementary experiment in which 
we adiabatically dress the rotating BEC with the RF field. 
In presence of the RF field of frequency u, the effective 
magnetic field is P e ff = (2^/^X0, f^R, —S), where (1b is 
the Bohr magneton and 6 = uo — iibB/(2%) is the detun- 
ing from resonance. On resonance, B e ^ oc y. At t = 
we adiabatically (in 100 ms) sweep 6 from a large value 
(^> Qr) to zero, thus preparing a P z = superposition state 
\V) = (It) + * ID) /V2 [Fig.gc)]. At this point, |?/) is equiv- 
alent to the \tt/2) state prepared by an RF pulse, which does 
not show long-term current stability [Fig.[2jb)]. However, if 
we leave the RF field on during the in-trap evolution, \y) is 
an eigenstate of the Hamiltonian and the coherence between 
|t) and \\) components does not decay [24]. In this case 
P z = supercurrent is stable and persists for more than a 
minute [Fig.gfd)]. 

These experiments clearly show that for analysing current 
stability in a partially polarised gas we must distinguish P z 



and |P|. With this understanding, we now quantitatively 
characterise the onset of the supercurrent decay in Fig. [3] by 
the time r at which the probability that the first phase slip 
(g = 3 ^ 2) has occurred is 50 %; this closely corresponds to 
the border of blue- and white- shaded regions. 

In Fig.^a) we see that r rapidly increases for P z > 0.64, 
saturating at 100 s due to the TV-decay. We now combine our 
measurements of r and the transverse-polarisation decay f(t) 
[Fig.^b)] to calculate |P| at the onset of supercurrent decay: 

P(r) = y/Pi + (1 - PJ)/(t) . (1) 

In Fig. [5jb) we clearly distinguish two regimes: one where 
P(t) is constant (within errors) and one where P(r) = P z . 
We thus complete our physical picture and accurately deter- 
mine the critical spin polarisation P c : 

(1) If P z > P c , then |P| never drops below P c , the super- 
current is fundamentally stable, r ^> t co h and P(r) = P z . 

(2) If P z < P c , supercurrent decay starts at r < t co h, when 
the decaying P becomes equal to P c . From all the data in this 
regime we get P c = 0.64(1). 

For < P z < P c the value of r varies from 4 s to 15 s 
and the orientation of P(r) = {y/P% — P|, 0, P z ) in spin 
space varies from purely transverse to purely longitudinal, but 
the onset of the supercurrent decay always occurs at the same 
|P| . We thus conclude that the region of supercurrent stability 
is in fact the outer shell of the Bloch sphere where |P| > P c 
[Fig.[5jc)]. This spin-rotational symmetry is intuitive but we 
note that it need not be universal. In our 87 Rb gas the strengths 
of intra- and inter-component interactions are almost identi- 
cal, so the Hamiltonian is almost invariant under rotations in 
spin space. In the future it would be very interesting to study 
supercurrent stability as a a function of both intra- and inter- 
component coupling strengths. 

The existence of a critical population imbalance for su- 
perflow stability was predicted in Refs. fT5lfT7lL assuming 
equal intra- and inter-component interactions and no inter- 
component coherence. The current instability was associ- 
ated with out-of-phase density fluctuations in the two com- 
ponents. However the agreement on the value of P c has not 
been reached. In Refs. (T5] [TSI it was predicted that any 
q > 1 flow is unstable for essentially any P < 1, but accord- 
ing to Ref. ifTTl such current is stable above some non-trivial 
interaction-dependent P c . The latter conclusion qualitatively 
agrees with our observations. However, none of the exist- 
ing theories is quantitatively applicable to our experiments, 
since they are limited to the simplified cases of reduced di- 
mensionality and very weak interactions. Moreover, the inter- 
play of the spin and rotational degrees of freedom may involve 
new physical effects. Specifically, the dynamics of the local 
spin vector on the Bloch sphere can result in a Berry phase 
and unwind the "scalar" phase describing the rotational flow; 
the timescale for this process would be the same as the spin- 
decoherence seen in a Ramsey experiment l25ll . 

In summary, we have observed persistent currents in multi- 
ply connected spinor condensates, demonstrated the existence 
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FIG. 5: (color online) Critical spin polarisation P c . (a) Characteristic time of the first phase slip, t, versus P z . Vertical blue line marks P c , 
accurately determined in (b). (b) \P\ at the onset of the supercurrent decay. Fit to the points with P(r) > P z (horizontal blue line) gives 
P c = 0.64(1). (c) Stability diagram on the Bloch sphere. The blue-shaded region in Fig.plmaps into the outer shell \P\ > P c . 



of a critical spin polarisation for stable superflow, and eluci- 
dated the role of spin coherence in supercurrent stability. Our 
results should stimulate further theoretical work on this fasci- 
nating many-body problem and are also relevant for applica- 
tions in trapped-atom interferometry. An important next step 
would be to study supercurrents in a two-species system with 
significantly different intra- and inter-component interactions. 
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